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Abstract. The stability of the orbital motion of two long cylindrical magnets 
interacting exclusively with magnetic forces is described. To carry out analytical 
studies a model of magnetically interacting symmetric tops [1] is used. The model was 
previously developed within the quasi-stationary approach for an electromagnetic fi- 
eld based on the general expression of the energy of interacting magnetic bodies [2]. 
A special role in the investigation of the stability of orbital motions is played by 
the so-called relative equilibria [3], i.e. the trajectories of the system dynamics whi- 
ch are at the same time one-parameter subgroups of the system invariance group. 
Nowadays their stability is normally investigated using two similar approaches - 
energy-momentum and energy-Casimir methods. The most suitable criterion for the 
system stability investigation was formulated in the theorem of [4]; this stability 
criterion successfully generalizes both the methods mentioned above and covers the 
Hamiltonian formalism based on Poisson structures [1]. The necessary and sufficient 
conditions for the circular orbit stability are derived from this theorem. 

1. Introduction 

The paper describes the new results on the research into the orbital motion stabi- 
lity of magnetically interacting bodies. It continues the cycle of papers devoted to the 
investigation of contact-free confinement of bodies by means of magnetic forces. 

There are two different historical aspects of our problem, physical and mathematical. 

The physical aspect is attributed to V.V. Kozoriz who was the first to come out in 
1974 with a supposition about the possibility of stable dynamic states in the so-called 
"compact" magnetic configurations. 

At that time an opinion about the global instability of electric and magnetic systems 
within the classical electrodynamics was widespread among physicists [5]. These ideas in 
statics were associated with Earnshaw's theorem; in dynamics they are associated with 
the so-called "l/r'^ problem". 

The supposition of V.V. Kozoriz seemed to be paradoxical and received no support in 
the scientific community; therefore he actually carried out his research on his own. 

In his papers, among which monograph |^ is worthy of mentioning about, he gave a 
whole number of models and predicted the orbital motion stability for them. 

Applying the physical idea about the nature of stability in such systems and analyzing 
the behaviour of the potential energy of the considered systems he gave a number of 
stability criteria for these systems. 

However, the method he used for investigating the stability cannot be called rigorous. 
Now it is quite clear that the proposed orbit stability conditions are not sufficient. 

Comparatively not long ago V.V. Kozoriz and his co-authors made numerical modelli- 
ng of the system of two long cylindrical magnets using a PC-cluster from the Institute 
of Cybernetics; they managed to discover a number of stable orbits [7J. However, they 
didn't produce any analytical proof of stability in this paper. 
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At the moment the hypothesis about the possibihty of stable orbital motions in 
"compact" magnetic systems was made the suitable methods of analytical stability investi- 
gation were only at the initial stage of their development and were unknown to physicists 
for a long period of time. 

This second, mathematical, aspect is associated with such names as B. Konstant, J.- 
M. Souriau, V. I. Arnold, A. A. Kirillov, and especially with J. E. Marsden. Marsden's 
lectures given in the Royal Mathematical Society and published in a separate book [3] 
are completely devoted to the stability problem. For several decades his scientific school 
has been developing the theory of stability of Hamiltonian systems with symmetry and 
applying it in the research into fluid flow, stability of plasma, elastic bodies, in the general 
theory of relativity and the quantum field theory. 

However, among the numerous results of these books and papers of his colleagues we 
could not discover any result which would be completely suitable for our system. This 
is due to the fact that the basic results on stability of Hamiltonian systems concern 
symplectic manifolds. The paper of J. -P. Ortega and T. S. Ratiu [4J is an exception; their 
theorem concerns an analytically more complicated case of Poisson structures. It is the 
theorem that we use in our paper. 

2. Hamiltonian formalism based on Poisson structures 

The Hamiltonian formalism based on Poisson structures (PS) [Ij gives an algebraic 
coordinate-free representation of dynamics which is especially important when a system 
consists of rigid bodies. 

In our case the Poisson manifold is the direct product of Euclidean spaces 



Thus, the generatrices of our dynamic system will be Xi as relative coordinates of 
two bodies; pi as pulse components of relative (orbital) motion; m^, n.i as components of 
angular momentums of the 1st and 2nd body, respectively; /Xi, Vi as components of unit 
vectors of symmetry axes for the 1st and 2nd body, respectively. 

The nonzero Poisson brackets between the generatrices on P are as follows, 



and the remaining influential Poisson brackets are equal to zero. 

It can be easily checked that the Casimir functions of this Poisson structure will be 



P = RlxRlxRlxRlxRlxRl 



(2.1) 



(2.2) 




{fl , m) = A/3 ~ consti, [v ,n) — N3 ~ const2. 



The Hamiltonian of the system can be written in the following form 



h ^ T + U(r,c ,c ,c ), 



(2.3) 



where 



r = |af|; e — x/r; 



< 



c = {e,fl); 
c" = (e, t7); 
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This form of the potential energy U describes a wide enough class of paired interactions 
of magnetic bodies such as rotationally symmetric permanent magnets, superconducting 
rings and solenoids, dipoles, and mixed systems [2j. 



3. Problem on stable orbital motion of two cylindrical 
magnets 

We will use the constructive proof method to prove the existence of stable orbital 
systems interacting with magnetic forces. It means that we will demonstrate the possibility 
of stability on the example of a specific physical system of this type. 

We will take two identical long cylinders as magnetic bodies. As it was shown earlier 
in our paper [1 : , the interaction of such bodies can be described by the Coulomb potenti- 
al energy of four fictitious magnetic charges [8J concentrated on the extensions of the 
cylinders. Thus, the potential energy of the system takes the following form 

y ee_ 

e',e"=±l ^ " 

where R^'^"{r,c ,c ,c ) = \/r^ + l'^ + l"^ + 2r{e"l"c" — e'l'c) — 2e' e" 1' I" c" ; I , / are 
semilengths of the cylinders; r is the distance between the centers of the cylinders; k , k 
are quantities of fictitious charges. 

The potential energy in point p.ip is a special case of the potential energy in (12.31) . 

Let us consider the problem in the center-of-mass system of the bodies [IJ. 

The kinetic energy of the system will be formed by the kinetic energies of translational 
(orbital) motion of the bodies and their self-rotation kinetic energies: 

rji/ 2 ^2 ^2\ 1 2 I '-'^ ^2 I ^2 

Tip,m,n)^—p +-m + -n , 

where r, p are orbital coordinates and pulses; m, n are angular momentums of the 1st 
and 2nd body, respectively; m — ^^^^ a ~ (3 — (if the principal moments of inertia 
of the body are equal to each other, Ii = I2 = I±, the rigid body is termed a symmetric 
top). 

Remark: If the consideration is limited to symmetric tops it is possible to use the 
uniform Cartesian frame, the components of physical vector quantities being suitable 
generatrices for the Poisson structure. In this case the Poisson brackets of dynamic vector 
variables can be considered, and the equations of motion can be written in a vector 
form ^Ij. 

The corresponding Hamiltonian equations of motion are as follows 
r = —p; 

^= -drUe- i(9,'[/Pl(M) +9,"C/Pl(z7)); 
/2 = a(jh X /i); 

m — d^'U{ex ft) - d^"'U{fI x v); 
a — I3[n X v); 

ft — d^" C/(e X z7) + d^'" U{fl x v), 

where operator Pj^ is the projection operator on the plane perpendicular to the vector e, 
i.e. Plijl) ^ fl — c e. 

It can be directly checked that the total momentum of j = ^ + to + n, I = x x p, for 
the Hamiltonian function (j2.3p is an integral of motion: j = {j, h} = 0. 
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4. Relative equilibria 



During the last decades large progress in the investigation of the system dynamics 
was reached thanks to the application of group-theoretic methods. This also concerns 
the research into the stability of some types of motions. In particular, a special role 
in the investigation of the stability of orbital motions is played by the so-called relati- 
ve equilibria i.e. the trajectories of the system dynamics which are simultaneously 
one-parameter subgroups of the system invariance group. Nowadays their stability is 
normally investigated using two similar approaches - energy-momentum and energy- 
Casimir methods. 

We will show that there are relative equilibria in the system in point. 

Let z be a fixed axis. Let us consider a subgroup of rotations around this axis. Each 
one-parameter subgroup of this group will be characterized by its angular velocity of 
rotation uj = loCz- The rate with which any physical quantity t; in our problem changes 
along the orbit of the given subgroup will be set by the formula v = Co x v. 

Therefore, the following relations should be fulfilled for the relative equilibrium 

|f=a;(e^xr); fl = uj{ez x fl); v ^ uj{ey^ x v)\ ^^^^ 
\v = ^{(^z 'X p)', m — uj{ez X rn); n — loiez x n). 

The fixed point on the orbit can be chosen as 




fl = 63; V = -63; 
m = mea; ft = nea, 



(4.2) 



where ei, 62, eij is a fixed basis. 

Thus, the bodies gyrate in the plane xy; the axes of the magnets are parallel to the 
axis z and are opposite in direction; and all moments, i.e. the orbital moment I — x x p, 
intrinsic moments m and n, and the total momentum j, are parallel to the axis z. 

So (|4.ip will then have the following form 

if=uj{ezxr); fl = 0; v = 0; ^^^^ 

This is in an agreement with the equations of motion at the given starting conditions 
since c — 0, c — 0, c =—1, and the evaluations show that in this case 9^'J7 = n 
d^"U = and 9^"C/ = 0. 

It means that the conditions of (|4.3p are reduced to the equality of centripetal and 
centrifugal forces. 

Remark: It should be noted that according to the conditions of the theorem of section 
5 it is required to prove a weaker statement, namely, the fulfilment of the conditions of 
in one point. 



5. Theorem of stability of relative periodic orbits in 
Hamiltonian systems with symmetry 

We will consider theorem 4.8 of paper [4J used to prove the orbital motion stability in 
our problem in more detail. We will quote its original statement: 

"Theorem 4.8 (Generalized energy- momentum method). Let (M, {,},G, J : M 
g* ,h : M ^ R) be a Poisson system with a symmetry given by the Lie group G acting 
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properly on M. Assume that the Hamiltonian h e C°°(M) is G-invariant and that J is 
equivariant. Let m G M be a relative equilibrium such that J{m) = fJ. d is compact, 

H := Gm, and ^ G Lie{NG^{H)) is its orthogonal velocity, relative to a given AdH- 
invariant splitting. If there is a set of G^-invariant conserved quantities Gi, C2, . . . , G„ G 
G°°(M) : M ^ R, for which 

- J« + G^ + G^ + • • • + G")(m) (5.1) 

and 

d^{h - J« + G^ + G^ + • • • + G")(m)|M/xM/ (5.2) 
is definite for some (and hence for any) subspace W such that 

kerdC\m) n • • • fl A:erdG"(TO) n fcerr„, J = 1^ © T,„(G^ • to), (5.3) 

then TO is a G^-stable relative equilibrium. If diniW = 0, then m is always a G^-stable 
relative equilibrium." 

Remark: The G^ - stability appearing in this theorem is given by the following defi- 
nition: 

"Definition 4.6. Let (Af, {, }, /i, G, J : M — g*) be a Hamiltonian system with 
symmetry and let G be a subgroup of G. A relative equilibrium to G M is called 
G -stable, or stable modulo G , if for any G -invariant open neighborhood V of the orbit 
G • TO, there is an open neighborhood [/ C y of to such that if Ft is the flow of the 
Hamiltonian vector field Xh and u ^ U, then Ft{u) G V for all t > 0." 

It is evident that there are five objects appearing in this theorem, i.e. M, {,}, G, J , h, 
where the two first objects M, {, } actually define the Poisson manifold (phase space), G is 
the Lie group acting on this manifold, J is the momentum map [9], h is the Hamiltonian 
of the system. 

The given theorem includes conditions of topological, algebraic and analytical character. 

The topological condition the "Lie group G acting properly on Af is of a technical 
nature and is always fulfilled for compact groups. The condition "G^ is compact" is 
essential. 

The algebraic conditions are requirements of invariance relative to an action of the 
group G: the Poisson structure, Hamiltonian function, momentum map (equivariance). 

Another algebraic condition defines the "orthogonal velocity, relative to a given Adn - 
invariant splitting". This notion is complicated enough and is the subject of the previous 
research of these authors. But everything is simplified in our case, and this quantity is an 
angular velocity parallel to a selected value of the momentum. 

To formulate the analytical conditions, a set of integrals of motion Gi,G2,...,G„ 
is introduced, and some efficiency function of the Hamiltonian and integrals of motion 
including moment components (in our case there is one component) is formed. Concerning 
Gi, G2, . . . , G„ there is another algebraic condition we did not mention before. This condi- 
tion is that these values should be invariant with respect to the subgroup G^. Equations 
(j5.m5.2p are similar to the relations for seeking a conditional extremum using the method 
of undetermined Lagrangian multipliers. 

However it is necessary to mention that, first of all, this similarity is not full, and, 
secondly, the theorem only requires checking conditions (j5.m5.2p in a certain point of 
orbit. 

Equation (|5.ip is analogous to the necessary condition for extremum, i.e. all partial 
derivatives of the efficiency function should be vanished in the point under test. 

Equation (j5.2p is more complicated and resembles the sufficient condition for extremum. 
We should establish positive definiteness of the efficiency function relative to only some 
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variation subspace W rather than to all possible variations of variables. Variations from 
W are characterised by the following properties: 

- they should conserve all integrals of motion involved in the efficiency function; 

- they should be transversal to the orbit direction. This is the main idea of requirement 
(ESI). 

Hamiltonian systems require the improvement of the notion of stability as for example 
the notion of asymptotic stability is not applicable to them. Indeed, let us suppose there 
is a stable orbit, as well as another orbit which is so close to the stable orbit that the 
stability is not broken, then for a rather long time interval the initially close points on 
these orbits can diverge quite considerably; they can for example be found on the opposite 
sides of the orbits. This behaviour is the case for Hamiltonian systems which could be 
well observed during the numerical modelling of our system [IJ. 

This well-known fact requires a suitable definition of the orbital stability which is given 
in paper [3] (see Definition 4.6 above). This definition contains the same system with 
symmetry as in the theorem. The stability is formulated with respect to some subgroup 
G . For this purpose a tubular neighbourhood of the stable orbit is introduced such that it 
entirely consists of the orbits of subgroup G . Then the trajectory of the system beginning 
in some neighbourhood of a stable orbit point should not go beyond the given tubular 
neighbourhood. 

6. Fulfilment of theorem conditions in the given problem 

First of all, it should be demonstrated that the group action is a Poisson action. 
The conceptual definition of a group action on the Poisson manifold and the associated 
formulas are given in book ^J. Below, the notations introduced in this book are used. 

$:GxP^P, $g(p)=.g-p. 

Actions <i> of the Lie group G on the Poisson manifold have the following form 

<^>;{f,,F2}^{^;f,,^;f2}. 

The action of the Lie algebra g of the Lie group G on the manifold P is then defined in 
terms of vector fields £,p{z) according to the formula: 

The application of modern group-theoretic methods in the investigation of our 
system is specific because these methods are far-reaching generalizations of the angular 
momentum theory whereas it is the angular momentum that is the momentum map in 
our case. 

Thus, in our case it is necessary to bring the complications and niceties required for 
the transition from a prime model to more and more general and complicated models 
down to the initial prime model of the rotation group actions and the associated angular 
momentum theory. 

Thus, it appears that our model is trivial neither mathematically nor physically. 

The group of transformations we are interested in is 5*0(3), i.e. a group of self-rotations 
of the Euclidean space, or a group of orthogonal 3x3 matrices with det = 1. This group 
is connected. 

The elements of the Lie algebra so(3) can then be considered either as antisymmetric 
matrices cj G so(3) or as vectors w e i?'^, i.e. iJik — eukoj' — —e-ikiLo^ — > ijJi = — ^^i^sWrs 
whence uj[v\ — uj x v, where v ^ . 
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The elements j G so(3)* of the space dual to so(3) can also be represented through 
three-dimensional vectors using the ordinary Euclidean scalar product as canonical pairing 
(J, uj) = J ■ Co. 

The coadjoint action of the group G = SO{i) in so(3)* will be reduced to the rotation 
in the Euclidean space 

Ad\^^=A^Ad\ = A-^, AeS'0(3). (6.1) 

In our case the Poisson manifold is a direct product of the Euclidean spaces (j2.1l) , and 
the Poisson brackets on P look like (j2.2p . 

The group action on P is reduced to a standard action of 50(3) on the Euclidean 
space, i.e. if a A G 5*0(3), and iT is a vector which is one of the 6 factors of (|2.ip . then 

HA,{.. V, ...)) = {.. .,A[vi...). (6.2) 

Physically it represents the rotation of the system as a whole; all physical quantities 
of a vector character are rotating synchronously. 

We will show, that the relations in (|2.2p are invariant relative to this action of 50(3) 
on P. For example, (57i,572 are constants), 

{■mi,mj] = Eijkrrik — > {D,\mi,nlmj} = Ql^l^^^.^jhrnk 

{(ill, to), {^2,m)} = (di X n2,rn). 
Similarly, the first relation in (j2.2p can be written as 

{(e,f),(/3,p)} = (e,^). 

It means that relations given below are equivalent to the corresponding base relations 
(|2.2p however they are invariant with regard to rotations: 

'{{^^x),0,p)} = {^J}; 
{{n,rn), {ff,fl)} ^{dx Tf,fl); 
< {{ni,m),{n2,Tn)} = {nixn2,m); (6.3) 
{(Il,n},(r7,z7)} = (fixrf,i7); 

where il, fii, i72, , t?, /? are some constant, i.e. independent of the generatrices, vectors. 

It follows that having transformed the base relations of (j2.2p into (j6.3p regardless of 
the coordinate system selected we actually proved that the action of the group 50(3) is 
a Poisson action on P. 

Now the field wp in our case will be written as 

wp(...,w, ...)) (...,w X w, ...). 

Let us show that the total angular momentum j in our problem is the momentum 
map. It means that it is necessary to specify the Hamiltonian for each field wp, i.e. 
Up = {z, Hi^} where z G P. 

Direct checking will demonstrate that the total momentum of the system j = x x p + 
rh + fi generates all such Hamiltonians H^^ = (w, j), i.e. the dynamic variable j is the 
total momentum of the system as well as the momentum map if so(3)* is equated with 
the Euclidean space R^. 
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The global equivariance generally means the following 



,-1 



o J. 



(6.4) 



Comparing the formulas in (j6.1l6.2l6.4|) we determine that in this special case the 
momentum map is also globally equivariant which is a rule for compact connected groups 



Therefore the following conditions are fulfilled in case of the dynamic system under 
consideration: 

1. The action of the group S0{3) and its Lie algebra so(3) on P is a Poisson action, 
i.e. it conserves the Poisson structure (Poisson brackets) on P. 

2. The momentum map is (globally) equivariant. 

3. The Hamiltonian of the problem and Casimir functions are invariant with regard 
to action (j6.2p of the group 5*0(3). 

4. The subgroup in terms of the theorem in section 5 is a one-parameter group of 
rotations around some fixed axis; this group is compact. 

5. The subgroup H — Gm in terms of the theorem in section 5 (Gz^ in our notation) 
is a trivial group composed of one element, i.e. group identity. 

6. The subgroup Nq (H) appearing in the theorem of section 5 is nothing else than 
G^. Accordingly, the algebra Lie{Noi^{H)) is one-dimensional and consists of rotation 
angular velocities around the same fixed axis. 

7. The integrals of motion used, namely the Casimir functions and component of the 
momentum along the selected axis, are invariant. 

All conditions of the theorem of section 5 relating to the actions of groups and Lie 
algebras on a Poisson manifold are fulfilled for our dynamic system. 

7. Selection of orbit point and integrals of motion 

The theorem of section 5 allows reducing the investigation of the relative equilibrium 
stability to the checking, in a fixed point of the orbit, of relations very similar in form 
to those occurring when checking the constrained minimum (maximum) by the method 
of undetermined Lagrangian multipliers. The necessary and sufficient conditions for the 
circular orbit stability are derived from this theorem. 

According to this approach the efficiency function is written as 



where w, A^, A'^, A'^, are Lagrangian multipliers, Ci are Casimir functions, and j is the 
momentum map corresponding to the Poisson action of the rotation group G = S'0(3) of 
the system as a whole. Accordingly, the equation for determining these constants is 



and the sufficient condition for the minimum is a positive definite quadratic form 



and, secondly, the subspace of variations 5z must be transversal with regard to the di- 
rection of motion along the orbit in the given point Zg- 



[9]. 



H = H-uj^ + A^Gi + A^Ga + A^Gg + X^Ci, 



dH, =0, 



(7.1) 
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The Casimir functions in our case are written as 



C4 — (iy, n). 



We select the point in a way shown in (|4.2I) . jo = j{ze) = xq x ^0 + m + n wiU then 
be directed along the axis 63. The stationary subgroup of this value of the momentum will 
be Gj^, i.e. the rotation subgroup around the axis 63. This group is compact as supposed 
in the theorem of section 5. Other premises of the theorem having a technical character 
are also fulfilled in our case (see section 6, items 1-7). 

10 variations of the following form can be considered independent 

Sx^,6x'^; S/i^ ,Sii^ ,6mi,Sm2', Siy^ ,6^^ ,Sni,Sn2- 

To meet the requirements mentioned above, the remaining variations are expressed in 
terms of the independent variations in the following way 





0; (5m3 = 0; Su^ = 


Sn3 = 


6x^ = 


{Smi + Sni)/po; 




< Spi = 






Sp2 = 


ro 




^5p3 = 


{Sm2 + Sn2)/rQ. 





0; 



(7.3) 



From (|7.1I) for Lagrangian multipliers w, we obtain 

Po/M = ujro; 
Pquj drll ; 
X2 — {uj ~ am)/ fi; 
A4 = (w — I3n)/v; 

Ai = - [vd^'" u + 111X2) /m; 

, A3 = - U + nA4) /v. 



If we substitute the expression for dependent variations (|7.2p into the quadratic form 
(Hessian) from condition (17. 3p we obtain a quadratic form from only 10 independent 
variations. Its positive definite will be the orbital stability as it is follows from the theorem 
in section 5. It appears that in the case of the following order of independent variations 

Sx^^Sx^] S^i^ ,5v^ ,5mi^5ni] Sfi^ ,6^^ ,Sm2,Sn2 

the matrix of this quadratic form takes a block-diagonal form, one 2x2 block and two 4x4 
blocks: 
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Thus, each of these blocks should define the positive-definite form from the correspondi- 
ng variations. The positive definite conditions are analytically deduced in the Maple 
system using Sylvester's criterion. Ten independent conditions for the system parameters 
were obtained which make it possible to define the regions of the orbital stability in the 
parameter space of the system. Most inequalities derived are quite complicated analytical 
expressions, and their complete analysis can only be conducted numerically. However, it 
is easy to show the system parameters for which these conditions will be automatically 
fulfilled. For example, the diameter and length of cylindrical magnets: d = 0.0025 [m], 
h = 0.02 [m]. The reduced mass of two rigid bodies m = 0.0003828816 [kg]. The values 
of the system parameters: a = 3.87228183489- lO"^ [kg-^m-^]; = 0.15546875 [A m^]. 
Intrinsic angular momentum of rigid bodies |m| = |ri|= 5 • 10^^ [m ^- kg- sec~^]. The rigid 
body orientation at the initial instant of time is characterized by conditions (j4.2l) . The 
radius of the stable orbit is then i?orb=0.01 [m], the impulse of force is Porb=0. 0006491 
[N-sec], and the revolution period is Torb=0. 037062129 [sec]. The stability conditions for 
the orbit were fulfilled as well as derived in the Maple system (see Fig.). 
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Fig. Test the fulfilment of the theorem conditions for specify orbit 



Summary 

The paper is devoted to the analytical investigation of the possible existence of orbital 
motions in the mathematical model [1] which we proposed earlier and which describes the 
magnetic interaction of rigid bodies. 

It has been shown that thanks to the application of group-theoretical methods the 
Hamiltonian formalism developed in [IJ on the basis of Poisson structures is a suitable 
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mathematical instrument for the analytical investigation of the orbital motion stability. 

The possibility of stable orbital motion has been proved for a specific system. 

Ten conditions for the system parameters have been analytically obtained; they make 
it possible to define the areas of the orbital stability in the parameter space of the system. 

Following the logic of the constructive proof it is possible to assert that magnetic 
systems can form stable orbits. This destroys the widespread opinion ^ about the global 
instability of magnetic systems. 
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